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Cygnus X-l: A Case for a Magnetic Accretion Disk? 


Michael A. Nowak 1 , B. A. Vaughan 2 , J. Dove 3 , and J. Wilms 4 

Abstract. With the advent of RXTE, which is capable of broad spec- 
tral coverage and fast timing, as well as other instruments which are 
increasingly being used in multi-wavelength campaigns (via both space- 
based and ground-based observations), we must demand more of our 
theoretical models. No current model mimics all facets of a system as 
complex as an x-ray binary. However, a modern theory should qualita- 
tively reproduce — or at the very least not fundamentally disagree with 
- all of Cygnus X-l’s most basic average properties: energy spectrum 
(viewed within a broader framework of black hole candidate spectral be- 
havior), power spectrum (PSD), and time delays and coherence between 
variability in different energy bands. Below we discuss each of these basic 
properties in turn, and we assess the health of one of the currently pop- 
ular theories: Comptonization of photons from a cold disk. We find that 
the data pose substantial challenges for this theory, as well as all other 
currently discussed models. 


1. Energy Spectra 

Historically, energy spectra have been the primary focus of both observers and 
theorists. The low state spectrum of Cyg X-l (see Fig. 1) for the most part is 
well-fit by a power-law with a photon index ~ 1.7. Additional facts to keep in 
mind are: the total luminosity of Cyg X-l is likely <0.1 f/Eddi that most 
other black holes show a strong soft component and softer power-law tail above 
> 0.1 L E dd (cf. Nowak 1995). The softer, high state of Cyg X-l is of roughly 
comparable luminosity to the low state (Zhang 1996). The soft components that 
do exist in the low state are weak compared to the power-law. 

The most popular models for this state are hot, optically thin models [with 
advection dominated disks (ADD) being one example, cf. Narayan et al. 1996], 
and Comptonization of seed photons from a cold disk (cf. Haardt & Maraschi 
1993; Dove, Wilms, & Begelman 1996). ADD theories provide an explanation 
for the hard to soft transition at ~ 0.1 LEdd, although they have yet to self- 
consistently calculate the (weak) observed features that Comptonization mod- 
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Figure 1. Best-fit corona model to unfolded, non-simultaneous Cyg 
X-l spectrum. Data are from BBXRT, TTM, (simultaneous with) 
HEXE, and OSSE (Dove, Wilms, & Begelman 1996). 


els attribute to reflection. One theoretical motivation for the Comptonization 
models is the recent work (cf. Balbus, et al. 1996; Hawley et al. 1996) demon- 
strating the existence of a powerful magnetic field-driven shear instability that 
is strongest in low-density regions. Comptonization models assume that most 
of the disk energy is dissipated in a corona, with roughly half being reprocessed 
by the disk into soft photons (cf. Haardt & Maraschi 1993). 

Fitting the spectra of Cyg X-l above 2 keV typcally requires r ea ~ 0.2 
and kT e ~ 150 keV. With roughly half of the luminosity being reprocessed by 
the disk, and a large fraction escaping the optically thin corona, we expect a 
large soft flux at ~ 200 eV. This is not observed (see Fig. 1 and Dove, Wilms, 
Begelman 1996, where more refined models will be presented). In addition to not 
offering an explanation for the low-high transition at ~ 0.1 Lsdd, this is a severe 
problem for the corona models. The models do, however, correctly account for 
the Fe-line and reflection features, with ~ half of the soft 200 eV excess being 
due to reflection (as opposed to intrinsic disk emission). ADD models, with a 
similar geometry and a strong source of hard photons, may also overproduce 
soft X-rays if they attempt to model the observed Fe-line with reflection off of 
the optically thick, geometrically thin outer disk. 


2. Phase/ Time Lags Between Soft and Hard Variability 

Cyg X-l shows strong X-ray variability (rms fluctuations ~ 40%), with timescales 
ranging from milliseconds to tens of seconds. This variability is usually measured 
by Fourier transforming the observed time series and then forming the power 
spectral density (PSD), i.e. the squared transform. One can also form the cross 
power density (CPD) by multiplying together the complex transforms from two 
time series. The Fourier frequency-dependent phase of the CPD is known as 
the phase delay, #<*(/), from which one can form the time delay, = ^/ 27 r/. 
Cygnus X-l shows a characteristic time lag of the hard photons behind the soft, 
with Td ~ 2 (// 0.01Hz) -2 / 2 3 s (cf. Miyamoto et al. 1992). This corresponds to 
a nearly constant phase delay of ~ 0.1 rad. 
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Figure 2. Effects of Compton coronae on intrinsic phase lags, as a 
function of input blackbody temperature and mean scattering path 
length (A es = 10 6 , 10 7 , and 10 8 cm). From Nowak & Vaughan (1996). 


The only currently proposed ADD model explanation for this delay invokes 
thermal disturbances propagating from the outer disk edge toward the center (cf. 
Manmoto et al. 1996). This model requires, however, a very slow propagation 
speed ~ 10 _4 c in order to be in quantitative aggreement with observations (cf. 
Vaughan & Nowak 1996). Corona models explain the delay as possibly arrising 
from the magnetic flares that are allegedly energizing the corona (cf. Nowak 
1994). If the flare temperature increases throughout its evolution, hard photons 
naturally lag soft photons. However, diffusion of photon scattering times within 
the corona will lead to a u smearing out” of this intrinsic time lag, depending 
upon the spectrum of the seed photons, the corona temperature, optical depth, 
and mean free scattering path (cf. Fig. 2; Miller 1995; Nowak & Vaughan 1996). 
Current observations rule out large scattering paths (~ 10 8 cm) with very soft 
(~ 150 eV) inputs. Upcoming RXTE observations (obtaining phase lags in 6-8 
energy channels out to ~ 250 Hz) will provide more stringent constraints. 


3. Coherence Between Soft and Hard Variability 

The coherence function is essentially the normalized average amplitude of the 
CPD (cf. Vaughan & Nowak 1996, Nowak & Vaughan 1996). It is a measure of 
the fraction of the two time series that can be interpreted as time-independent 
and linear transforms of one another. Coherence must be less than or equal to 
unity, and be equal to unity only when there is a linear, time-indpendent transfer 
function relating one channel to the other. Coherence is less than unity under 
any of the following conditions: variability is due to thermal flares that are dom- 
inated by emission from the Wien tale; there are multiple, independent sources 
of variability (whether linear or not) with multiple, independent responses to the 
sources (whether linear of not); there is any non-linear transfer function from 
one channel to the other (cf. Vaughan & Nowak 1996). 

If coronae are truly formed via multiple magnetic flares, then we expect a 
strong loss of coherence, unless the the response to each of the flares is identical. 
Nowak (1994) showed that a kinematic model, with multiple flaring regions, 
reproduced both the PSD and phase lags seen in the very high state of GX339- 

4, as well as quantitatively agreed with the energetics of coronal formation. This 
model, and all others like it, produces a strong loss of coherence (Fig. 3). This 
is not seen in Cyg X-l (Fig. 3; Vaughan & Nowak 1996), which shows near 
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Figure 3. Left: Coherence of theoretical model with multiple flares 
(based on Nowak 1994). Right: Ginga data of Cyg X-l (cf. Vaughan 
& Nowak 1996) The Cyg X-l loss of coherence is likely a noise artifact. 


unity coherence between all well-observed energies. This also disagrees with the 
model of Manmoto et al. (1996). For that case the disturbance shocks upon 
reflection off of the disk inner edge, which is an inherently non-linear process, 
and thus leads to less than unity coherence for a stochastic source of variability. 

As they currently stand, coronal models of Cyg X-l have several problems. 
They predict too much soft X-ray flux, and they have yet to offer a viable 
explanation for the phase lags that also preserves unity coherence. ADD models 
also do not offer an explanation for the observed coherence, and they have yet 
to self-consistently calculate the observed “reflection” and line features. 
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